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I  SUMMARY  AND  CONCLUSIONS 

The  Ballistic  Missile  Defense  (BMD)  non-nuclear  kill 
problem  contains  impact  conditions  involving  relative  velo¬ 
cities  ranging  from  40,000fps  to  less  than  10,000  fps,  many 
possible  impact  angles  depending  on  end  game  geometries  and 
possible  multiple  impacts.  These  engagements  conditions  can¬ 
not  be  simulated  in  existing  experimental  facilities,  there¬ 
fore  a  credible  three-dimensional  (3D)  analytical  tool  which 
can  evaluate  the  impact  process  at  every  stage  through  its  time 
history  is  essential  to  better  understanding  the  impact  pheno¬ 
mena  and  resolving  the  non-nuclear  lethality  problem. 

No  single  numerical  code  exists  with  the  capability  of 
analyzing  the  complete  hypervelocity  impact,  the  penetration 
process  and  late  time  effects  of  material  failure  unique  to 
BMD  engagements . 

This  study  concludes  that  a  reasonable  approach  to  the 
analytical  problem  may  be  resolved  by  coupling  (output  from 
one  numerical  technique  becomes  the  input  to  another  numerical 
technique)  existing  numerical  techniques. 

Three-dimensional  Eulerian  hydrodynamic  codes  (coupled 
with  strength  of  materials)  are  available  for  the  early  times 
where  extremely  high  pressures,  non-linearities,  and  large  de¬ 
formations  exist.  Their  shortcomings  are  at  times  where  the 
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stresses  are  small  and  where  the  details  of  material  interfaces 
are  important.  On  the  other  end  of  the  deformation  spectrum, 
three-dimensional  finite  element  codes  capable  of  including 
non-linear  effects,  detailed  material  interface  information, 
and  small  deformation  with  reasonable  failure  models  are  also 
available. 

To  couple  techniques  such  as  these  in  a  dynamic  model, 

is  a  lengthy  (if  at  all  possible)  task,  that  would  require 
computers  of  the  size  not  yet  available  and  enormous  computer 

times. 

The  possibility  of  coupling( static)  these  techniques 
throughout  phases  of  the  time  history  of  the  impact  processes, 
seem  reasonable  and  achievable  in  a  short  period  of  time.  Some 
thoughts  on  coupling  are  presented  in  Section  IV. 

A  re-run  (Section  II)  of  the  NRL  45°  oblique  impact 
(ref  l)  has  been  completed  out  to  a  time  of  'V'  11.0^  sec. 
where  we  studied  the  effect  of  zoning.  Zoning  apparently  had 
little  effect  in  changing  the  peak  pressure  for  gauge  number  (l), 
the  closest  gauge  to  the  front  surface  of  the  target.  The  area 
under  the  curve  of  pressure  versus  time  appeared  to  remain 
constant. 

A  two-dimensional  calculation  of  NRL  #5  was  performed 
with  the  SOIL  code  for  the  purpose  of  using  output  from  the 
SOIL  code  as  input  to  a  structure  code.  This  calculation  and 
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results  are  reported  in  Section  III. 

Some  thoughts  of  adding  techniques  that  may  add  to  the 
accuracy  of  2  and  JQ  calculations  and  allow  larger  grids  are 
presented  in  Section  V. 

Some  ideas  as  how  to  begin  calculations  in  2  and  3  dimen¬ 
sions  and  accuracy  or  credibility  of  results  are  presented  in 
Section  VI. 
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II  RE-RUN  OF  45°  NRL  TEST  PROBLEM 

A  TRIDORF  calculation  of  the  45°  oblique  NRL  experiment 
was  reported  in  ref  (l). 

To  summarize  briefly,  the  peak  pressure  in  the  gauge 
nearest  to  the  front  surface  of  the  aluminum  block  (target) 

indicated  a*  52  kilobars  from  the  experiment  and  32  kilobars 

from  the  code  run. 

We  felt  that  there  was  some  uncertainty  in  the  calcula¬ 
tion  due  to  the  coarse  zones  used  in  the  3-dimensional  calcu¬ 
lation.  Due  to  this  uncertainty  and  also  the  need  for  early 
time  data  for  a  coupling  experiment,  v<e  decided  to  re-run  the 
calculation  with  finer  zoning. 

Actually  3  calculations  will  be  reported.  The  first  one 
is  the  original  calculation,  the  second  one  is  where  finer 
zoning  was  placed  between  the  OTWR  plate  and  the  aluminum 
block  and  for  the  third  calculation,  we  used  the  zoning  of  the 
second  but  zoned  the  aluminum  block  finer. 

Figure  1  indicates  the  zoning  in  the  original  and  re¬ 
run  problems.  Figure  2  shows  early  time  comparisons. 

The  pressure  versus  time  for  the  3  calculations  for 
gauge  number  (l)  are  displayed  in  Figure  3. 

As  indicated,  the  finer  zoning  between  the  OTWR  plate 
and  target  resulted  in  a  slight  increase  in  peak  pressure  for 
gauge  no.  (l).  Using  the  same  zoning  between  the  OTWR  plate 
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and  target,  but  using  finer  zoning  in  the  target,  resulted  in 
an  additional  slight  increase  in  peak  pressure  for  gauge  no.(l). 

We  believe  that  any  additional  fine  zoning  will  not 
result  in  a  further  increase  in  peak  pressure  in  gauge  number 
(l).  These  last  2  problems  were  only  run  to  a  time  of  11. 

/u.  sec.,  due  to  the  limitations  of  core  size  available  for  the 
ARC  computer. 

A  zone  sensitivity  of  the  original  projectile  and  OTWR 
plate  was  not  initiated.  Even  though  the  computer  run  gives 
/\f  the  correct  degraded  projectile  velocity  and  projectile 
shape  (Figures  2  and  3)»  the  possible  artificial  diffusion  of 
the  projectile  (due  to  zoning  restrictions)  just  before  impact¬ 
ing  the  aluminum  target  may  be  responsible  for  the  code  pre¬ 
dicting  the  peak  pressure  for  gauge  number  (l). 

The  peak  pressures  for  gauge  number  (2)  and  (3)  agree 
very  well  with  experiment.  Another  possibility  is  that  the 
experiment  may  have  some  difficulties  for  gauge  number  (l). 

As  indicated  in  Figure  3  the  area  under  the  curve  (J  Pdt) 
is  fairly  constant  for  all  3  calculations.  Perhaps,  for  load¬ 
ing  of  inner  structures,  that  the  impulse  (integral)  is  only 
important  and  that  the  magnitude  of  the  pressure  pulse  is  not 
all  that  important.  Additional  discussion  of  this  idea  is  ela¬ 
borated  in  Section  VI. 
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ly  run  3D  calculation  can  give  reasonable  and  us 
to  later  time  structural  techniques. 

We  will  continue  these  discussions  in  Section  VI. 
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Ill  TWO-DIMENSIONAL  CALCULATION  OF  NRL  #5 
(CANDIDATE  FOR  COUPLING  EXPERIMENT) 

In  Section  II,  we  mentioned  the  possibility  of  using 
output  from  the  re-run  of  the  3D  calculation  as  input  for  a  3d 
finite  element  code.  Due  to  the  inability  of  acquiring  3D  fin¬ 
ite  element  codes  and  the  personnel  to  assist  in  the  coupling, 
we  chose  to  run  a  2D  calculation  and  consider  coupling  these 
results  with  the  EPIC  family  of  codes  (Gordon  Johnson  at  Minn¬ 
eapolis  Honeywell). 

Figure  4  displays  the  experimental  configuration  as  used 
in  the  code  calculation  (SOIL) . 

The  degradation  of  the  steel  projectile  passing  through 
the  heat  shield  material  agrees  very  well  with  the  experiment. 

The  code  -  experimental  comparison  data  is  displayed  in  Figure  4. 
As  can  be  seen,  the  calculated  pressures  are  consistently  lower 
(a  factor  of  2)  than  the  experimental  results. 

This  may  be  due  to  the  coarse  zoning  used  through-out 
the  calculation.  The  possibility  of  having  the  wrong  configura¬ 
tion  (density  and  shape)  due  to  zoning  between  the  plate  and  tar¬ 
get  could  also  explain  the  differences  observed. 

Massless  tracer  particles  were  placed  in  an  orthogonal 
grid  through-out  the  target.  The  location  of  the  tracer  parti¬ 
cles  (the  2  coordinates ) ,  the  pressure  energy  and  velocities  as 
a  function  of  time,  are  available  as  input  to  a  structure  code. 
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|  Again,  we  address  the  question  of  whether  impulse  is 

sufficient  or  whether  the  magnitude  of  the  pulse  is  required 
in  order  to  give  meaningful  input  data  to  a  structure  code. 
Again,  discussion  on  the  alternative  (fine  zones  with  an  in¬ 
crease  in  computer  time)  is  left  for  Section  VI, 

Vie  will  continue  to  use  this  2D  calculation  to  study  the 

,  effects  of  zone  size  and  material  properties. 

« 
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IV  COUPLING  IDEAS 


We  have  chosen  the  two-dimensional  world  as  a  beginning 
for  some  coupling  experiments.  The  results  (Section  III)  from 
the  SOIL  (2Dcalculation)  in  the  form  of  velocities  as  a  function 
of  time  and  2  coordinates  will  be  imposed  as  input  to  a  finite 
element  code.  The  objective  will  be  to  determine  if  the  velo¬ 
city  and  stress  field,  some  distances  into  the  target  will  dup¬ 
licate  the  SOIL  results  given  velocities  as  a  function  of  time 
at  a  moving  lagrangian  set  of  points  (close  to  the  initial  front 
surface  of  target).  This  coupling  experiment  is  presently  being 
worked  on. 

Other  potential  possibilities  may  be  to  apply  stresses 
at  this  moving  lagrangian  boundary  to  correlate  stresses  and 
velocity  fields  some  distance  into  the  target. 

We  feel  that  the  techniques  developed  for  the  2D  geometry 
will  be  applicable  to  the  JD  world  also. 

Failure  models,  fracture  and  spall  events  are  considerably 
easier  to  adopt  to  a  Lagrangian  or  Finite  element  technique  due 
to  the  ease  of  interface  (free  surface  or  material)  treatment 
in  techniques  where  the  mass  in  a  given  zone  is  preserved . 'These 
models  exist  in  rather  primitive  or  simple  formulation  as  of  now. 
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V  NEW  CONCEPTS  FOR  2D  AND  3D  TECHNIQUES 

Most  Eulerian  techniques  are  adopted  to  rectangular  zones 
rather  than  quadrilateral  zones.  As  indicated  in  Figure  5, 
quadrilateral  zoning  of  Shells  in  2  dimensions  or  generalized 
zoning  in  3  dimensions  may  be  sufficiently  more  accurate  and  re¬ 
quiring  fewer  computational  cells.  A  moderate  effort  has  been 
initiated  and  will  continue  during  this  year.  It  is  anticipated 
that  one  will  use  a  well  formulated  Lagrangian  code  for  the  hy¬ 
drodynamics  and  add  an  advection  model  with  the  appropriate 
geometric  modifications  for  fluxing  mass,  momentum  and  energy 
across  non-normal  interfaces  or  planes. 

Another  approach  enabling  one  to  have  a  larger  number  of 
computational  zones  (actually,  the  limitation  will  be  the  com¬ 
puter  time)  to  represent  the  problem  of  interest  is  the  pack¬ 
aging  concept.  Figure  5  is  a  possible  representation  of  the 
packaging  concept  (dividing  a  large  problem  into  N  smaller  prob¬ 
lems).  These  N  problems  would  communicate  across  a  common 
boundary,  allowing  the  zoning  to  be  different  in  each  of  the  N 
sub-problems.  This  would  allow  one  to  have  fine  zoning  only  in 
the  problem  that  requires  it,  rather  than  thru  the  entire  N 
problems,  A  version  of  a  2D  Eulerian  (RADOIL)  was  successfully 
modified  for  the  packaging  concept.  We  would  incorporate  as 
many  of  those  features  that  would  be  applicable  into  existing 
2D  and  3D  Eulerian  codes. 
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VI  SOME  IDEAS  AS  TO  ACCURACY  AND  ZONING 
FOR  2D  AND  3D  PROBLEMS 

As  indicated  in  Section  II  (3D  RE-RUN)  the  varying  of 
zone  sizes  did  not  effect  the  peak  pressure  for  the  gauge 
nearest  the  front  surface  of  the  target. 

However,  in  Section  III  (the  2D  calculation),  the  peak, 
pressures  as  calculated  were  consistently  lower  by  a  factor  of 
2  when  compared  to  the  experimental  results.  However,  the  vel¬ 
ocity  degredation  and  hole  size  in  the  target  correlated  very 
well  with  experimental  results.  Again,  we  postulated  that  finer 
zoning  might  bring  the  comparison  into  much  better  agreement. 

One  question  that  we  might  ask  is  would  this  change  the  impulse 
delivered  to  the  target. 

Zoning,  indeed  is  a  crucial  question.  A-prior,  do  I  know 
how  to  zone  up  a  2  or  3D  problem  such  as  to  remove  the  possible 
uncertainties  associated  with  a  finely  or  coarsely  zoned  prob¬ 
lem. 

Perhaps  so  systematic  approach  to  this  question  is  possible. 

l)  Run,  whenever  possible,  a  one-dimensional  calculation, 
approximating  the  pulse  to  find  a  consistent  set  of  zone  sizes 
that  would  produce  the  same  pressure  magnitude  and  pulse  shape. 
For  a  one-dimensional  slab  geometry,  the  only  attentuation  of  the 
pulse  will  result  from  the  rarefaction  from  the  bottom  surface, 
however  a  sperical  or  cylindrical  pulse  with  a  shell  of 
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approximate  thickness  in  the  2  or  3  dimensional  configuration 
would  produce  divergence  such  as  to  attentuate  the  pulse. 

2)  It  has  been  our  experience,  that  once  the  calculation 
(2  or  3D)  has  been  completed,  it  requires  an  extensive  period  of 
time  for  the  analysis  of  the  results.  It  would  be  more  practi¬ 
cal,  if  one  can  eliminate  some  of  the  uncertainty  in  results 
due  to  zoning,  before  going  ahead  with  the  large  full  calcula¬ 
tion. 

3)  Pine  zoning  may  be  required  at  spall  surfaces,  since 
spalling  is  a  function  of  the  peak  stress,  but  for  momentum  de¬ 
livered  to  a  target,  it  may  only  require  that  the  impulse  be 
calculated  properly. 

4-)  Any  discontinuity,  free  surface,  material  interface  or 
shock  front  will  require  a  minimum  number  of  zones.  For  example, 
3  zones  are  required  to  represent  a  shock  in  almost  any  numeri¬ 
cal  technique  and  approximately  the  same  number  of  zones  (mini-) 
mum)  will  be  required  to  describe  a  shell  or  layer  of  a  given 
material . 

5)  With  these  restrictions,  it  is  very  important  to  know 
whether  the  impulse  is  adequate  or  whether  the  pressure  (shape 
of  impulse)  magnitude  is  necessary. 

6)  For  a  2D  calculation,  requiring  zoning  that  is  a  fac¬ 
tor  of  2  finer  in  both  directions,  transposes  into  a  factor  of  8 
in  additional  computer  time. 
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For  a  3D  calculation,  this  could  result  in  a  factor  of  16 
in  computer  time.  It  may  be  that  one  would  calculate  the  same 
result  with  fewer  zones,  this  information,  at  the  beginning  of 
the  calculation,  would  be  extremely  useful. 

7)  The  configuration,  and  the  size  of  the  computer,  and 
time  and  funds  available,  may  dictate  the  zoning. 

An  example  of  this  approach  has  been  reported  in  a  three- 
dimensional  calculation  of  the  venting  process  for  an  -  on  axis 
explosion  in  the  MX  trench  (ref  2).  The  comparison  between  the 
finely  zoned  2D  and  the  3D  zoned  problem  (using  all  available 
core  of  a  CDC  7600)  is  remarkable.  The  one  dimensional  calcu¬ 
lations  would  still  give  one  some  feelings  to  the  credibility 
of  the  results.  These  results  are  shown  in  Figures  6,7  and  8. 

8)  An  impact  calculation  (ref  3)  of  a  cadmium  plate 
backed  up  by  a  void  and  a  iron  plate  gave  excellent  comparisons 
of  pressure  on  the  iron  plate,  hole  size  in  cadmium  plate  with 
experimental  results  (Figures  9  and  10).  The  final  calculation 
was  completed  only  after  a  series  of  preliminary  2D  calculations 
were  completed  investigating  the  pressure  pulse  in  the  first 
target,  and  the  zoning  required  to  adequately  represent  the 
debris  cloud  impacting  the  iron  plate.  The  number  of  zones  in 
each  plate  and  the  separation  between  plates  were  varied  in 
addition  to  the  material  properties. 
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FLASH  X-RAY  AT  1=4.5^  sec. 


Figure  2 

PROJECTILE  CONFIGURATION  as  a  FUNCTION  of  TIME 
and  an  EXPERIMENTAL  COMPARISON  at  t=4.5  sec. 
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Figure  4 

2D  CALCULATION  COMPARISONS  WITH  EXPERIMENT 
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PEAK  PRESSURE  ON  TUNNEL  AXIS  AS  A 
FUNCTION  OF  TIME  COMPARING  RESULTS  OF 
3-DIMENSIONAL  CALCULATIONS  WITH  THOSE 
OF  TRANSVERSE  AND  LONGITUDINAL 
TWO-DIMENSIONAL  CALCULATIONS 


Figure  7 

3D  and  2D  CALCULATIONS  of  PEAK  PRESSURE 
ALONG  TUNNEL  AXIS 
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CADMIUM  SPHERE  (R=.363  cm.) 
VELOCITY (7. 54  km/sec.) 


Figure  9 

CODE  CONFIGURATIONS  for  the  2  CALCULATIONS 
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